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Section  I. 


Introduction 


JMUttlrtlM. 


Overview 


Current  D062  safety  levels  computations  utilize  the  Laplace 
Distribution  to  approximate  the  distribution  of  demand  in  a  lead 
time.  This  paper  presents  several  plots  which  compare  the  Laplace 
distribution  with  an  empirical  model  of  D062  demand  in  a  lead 
time.  Section  I  provides  additional  background  for  the  paper, 
while  Section  II  presents  a  detailed  comparison  of  several  spe¬ 
cific  lead  time  demand  models. 


Background 

Safety  level  computations  utilized  in  the  Economic  Order 
Quantity  (EOQ)  Buy  Computation  System  (D062)  are  based  on  formu¬ 
las  originally  developed  by  Presutti  and  Trepp  (1970).  These 
authors  consider  the  problem  of  determining  order  quantities 
and  reorder  points  for  each  item  in  a  single-echelon,  multi -item, 
continuous  review  inventory  system  so  as  to  minimize  total  system 
holding  and  shortage  costs.  In  addition,  they  assume  there  is 
a  constraint  on  either  total  units  back  ordered  or  on  the  average 
number  of  units  in  a  back  order  position.  Presutti  and  Trepp 
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uted.  Howeve^,  jthey  then  utilize  the  Leplece  distribution  to 
jnonial.  With  this  subsitution,  Presutti  and 
Tvd^fP*  oUttfin  clqsed  fora  expressions  for  the  optimum  order  quantity 
and  reorder  point .  For  convenience,  we  refer  to  these  resulting 
fojraulas  as  the )PT> Formulas. 

Subsequent  simulation  studies  using  actual  demand  history 
for  Air  Force  items  showed  the  PT- formulas  were  significantly 
more  cost  effective  than  the  inventory  level  computations  then 
in  use;  that  is,  the  PT-formulas  provided  lower  levels  of  back* 
orders  for  a  given  investment  in  inventory  than  the  previous 
formulas,  or  conversely,  a  given  back  order  level  could  be  achieved 
with  the  PT* formulas  with  a  smaller  investment  in  safety  stocks. 

As  a  result  of  these  studies,  the  Air  Force,  the  Defense  Supply 
Agency,  and  the  US  Army  (for  high  demand  items)  adopted  the  PT- 
formulas  for  the  management  of  EOQ-type  items. 

Three  of  the  major  assumptions  which  are  embedded  in  the 
Air  Force  D062  implementation  of  the  PT-formulas  are  the  following: 

1.  Demand  in  a  lead  time  is  normally  distributed. 

2.  The  lead  time  is  known  and  constant. 

3.  The  mean  and  standard  deviation  of  lead  tine  demand 
may  be  accurately  estimated  from  available  history. 
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The  ebove  assumptions  are  commonly  employed  in  many  commercial 
inventory  systems,  and,  as  noted  above,  simulation  studies  have 
shown  the  resulting  formulas  are  significantly  more  cost  effec¬ 
tive  for  the  control  of  Air  Force  EOQ  inventories  than  the  pre¬ 
viously  used  formulas.  However,  several  recent  studies  have 
indicated  that  the  above  assumptions  may  be  a  poor  approximation 
to  the  actual  character ics  of  many  Air  Force  EOQ  items.  In  partic 
ular,  in  Reference  2  it  is  observed  that  the  distribution  of 
forecast  errors  appears  better  described  by  a  nonsymmetrical 
exponential  distribution  than  by  the  normal  distribution.  In 
Reference  4,  Hayya  observes  that  the  replenishment  lead  time 
for  many  D062  items  appears  to  be  highly  variable,  and  that  the 
limited  amount  of  data  on  procurement  lead  times  makes  accurate 
estimations  of  lead  time  parameters  difficult. 

As  a  result  of  these  findings,  we  have  used  historical  D062 
data  to  develop  a  refined  model  for  the  distribution  of  demand 
in  a  lead  time.  In  this  paper,  we  provide  comparisons  of  this 
empirically  derived  model  with  the  Laplace  distribution.  First, 
however,  let  us  consider  each  of  these  formulas  in  more  detail. 


The  Laplace  Distribution 

Let 


x  -  number  of  units  observed  during  a  procurement  lead  time 


I-4» 


* 
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u  •  the  expected  demand  in  a  lead  time. 

9  *  the  standard  deviation  o£  demand  in  a  lead  time. 

t  •  the  lead  time  in  months. 

With  the  above  definitions,  the  probability  density  function 
f(x)  for  the  Laplace  distribution  is  given  by 

Cl)  f(x)  -  _1 _  exp  (  -/Tiki  ) 

/T  or 

where 

(2)  k  ■  Cx-u)  /  <r 

i.e.  k  denotes  the  number  of  standard  deviations  that  the  demand 
value  x  exceeds  the  expected  demand  in  a  lead  time.  Given  (1), 
Presutti  and  Trepp  show  that  the  probability  that  the  demands 
actually  observed  in  a  procurement  lead  time  (X)  is  less  than 
or  equal  to  a  specific  numerical  value  x  is  given  by 


(3)  PCX  *  x)  • 


.5  exp  C  /T  k  )  for  k  £  0 

1.  -  .5  exp  C  -  /I  k  )  for  k  *  0 
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In  establishing  cost  effective  safety  stocks,  the  emulative 
distribution  function  PCX  *  x  )  is  particularly  important.  A 
common  approach  for  establishing  safety  stocks  is  to  consider 
the  trade-off  between  holding  costs  and  shortages  to  determine 
an  optimal  fill  probability  P*.  The  cumulative  distribution 
function  P(X  i  x)  is  then  used  to  determine  the  specific  value 
of  x  that  corresponds  to  this  optimum  fill  probability. 

As  noted  above,  an  important  assumption  embedded  in  the  D062 
safety  level  computation  is  that  the  parameters  u  and  <r  of  the 
lead  time  distribution  may  be  accurately  estimated  from  available 
data.  Since  demand  in  a  lead  time  is  not  directly  observed  in 
the  0062  system,  these  parameters  must  be  estimated  from  other 
data.  At  present,  the  following  estimation  equations  are  used. 
First,  let  R  denote  the  average  quarterly  demand  rate  observed 
over  the  past  eight  quarters,  and  let  QMAD  denote  the  Mean  Abso¬ 
lute  Deviation  CMAD)  associated  with  this  quarterly  demand  rate 
estimate.  Then  the  parameters  u,  <T  of  the  lead  time  demand  distri¬ 
bution  are  estimated  as  follows: 


u  «  R*t 
C4) 

&  -  0.5945  *  QMAD*  (  0.8235  +  0.42625  *  t) 

Where  ”*"  denotes  multiplication.  The  first  equation  is  derived 
from  the  fact  that  the  expected  demand  in  a  lead  time  of  t  periods 
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is  equal  to  t  times  the  expected  demands  in  a  single  period. 

The  standard  deviation  estimate  $  is  based  on  an  approximation 
suggested  by  Brown  (1967).  This  approxiaation  accounts  for 
the  fact  that  demand  rate  estimates  are  based  upon  averages  of 
random  variables  and  are  thus  correlated  from  period  to  period. 

An  Empirical  Model  of  Forecast  Errors  in  a  Given  Time 


Reference  2  presents  the  results  of  statistical  studies  to 
identify  the  actual  distribution  of  forecasting  errors  associated 
with  current  D067  forecasting  methods.  In  this  reference ,  actual 
CY71-79  D062  demand  histories  for  Sacramento  and  Oklahoma  City 
Air  Logistics  Centers  are  used  to  (a)  forecast  demands  in  a  given 
lead  time  using  current  D062  forecast  rules,  and  (b)  to  compute 
the  distribution  of  forecast  errors  associated  with  these  fore¬ 
casts.  Analytical  approximations  to  the  empirical  data  are  then 
developed.  As  a  result,  it  was  found  that  for  items  with  demand 
rates  of  three  units  per  quarter  or  more,  the  cumulative  distri¬ 
bution  function  for  demand  in  a  fixed  lead  time  of  t  periods 
may  be  approximated  by 

0.669  exp  (  0.7979  z)  for  z  *  0 

(5)  P(x  <  x  |  t)-  - 

1.  -  0.331  exp  (  -  0.463  z)  for  z  *  0 

where 
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(6)  t  -  (x  -  R  t  )/  (QMAD  /F  ) 

lliis  model  appears  to  be  a  particularly  good  fit  to  the  actual 
distribution  of  forecast  errors  for  positive  values  of  z. 


A  Bapirical  Model  with  Gama  Lead  Times 

Equation  (5)  describes  a  useful  model  for  the  distribution 
of  demand  in  a  fixed  lead  time  of  t  periods.  If  lead  time  is 
random  and  independent  of  demands  per  period,  the  unconditional 
distribution  of  demand  in  a  lead  time  may  be  found  by  averaging 
the  conditional  distribution  (S)  with  the  probability  distribu¬ 
tion  for  lead  time.  Specifically,  let  g(t)  denote  the  probabil¬ 
ity  density  function  of  lead  time.  Then  it  may  be  shown  that 
unconditional  distribution  of  demand  in  a  lead  time  is  given 


(7)  P(X  -  x)  - 


PCX  «  xlt  )  gCt)  dt 


In  Reference  4,  Hayya  describes  statistical  studies  to  iden¬ 
tify  an  appropriate  model  for  the  distribution  of  lead  times 
for  D062  items.  He  observes  that  several  probability  distribu¬ 
tions,  including  the  normal,  gamma,  exponential,  Weibull  and 
log  normal,  are  consistent  with  the  available  lead  time  data 
for  a  number  of  D062  items.  In  reviewing  Hayya' s  results,  the 
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Gama  distribution,  in  particular,  appears  to  be  a  useful  descrip* 
tion  of  lead  tine  for  the  purposes  of  this  study.  Specifically, 
if  lead  times  are  gaana  distributed,  we  have 

(8)  g(t)  -  J. _  ta_1  exp  (  -t  /  b  ) 

V(a)  ba 


where  a  and  b  are  parameters  of  the  distribution.  The  expected 
value  and  variance  of  the  Gamma  distribution  are  given  by 


(9)  ECt)  -  a  b 

(10)  Var(t)  -  a  b2  -  b  E(t) 


Hence,  one  method  of  establishing  the  parameters  for  a  Gamma 
distribution  is  to  estimate  the  mean  and  variance  of  lead  time 
from  historical  data,  and  then  use  (9)  to  solve  for  the  specific 
(a,b)  values  which  yield  the  desired  moments. 

A  second  estimation  procedure  is  based  upon  the  fact  that 
the  coefficient  of  variation  c  for  the  Gamma  distribution  is 
given  by 
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(11)  c  •  /Var(t)  ”  |/ab^  “  1 

E(t)  ab  /a" 

Hence, 

(12)  a  -  1_ 

c2 

Once  .a  is  known,  parameter  b  may  be  obtained  using  (9).  Specifically, 

(13)  b  -  E(t)/a 

Thus,  to  obtain  (a,b)  estimates  we  may  estimate  the  coefficient 
of  variation  c  and  the  expected  lead  time  E(t),  and  then  use 
these  values  in  (12)  and  (13)  to  obtain  (a,b)  estimates.  In 
the  calculations  reported  in  Section  II,  we  have  used  this  second 
approach. 


Estimates  for  a  Specific  Lead  Time  Distribution 

In  Reference  4,  Appendix  C,  Hayya  presents  estimates  of  the 
mean  and  coefficient  of  variation  associated  with  historical 
lead  time  data  for  62  EOQ  items.  Table  1-1  presents  a  summary 
of  the  coefficients  of  variation  observed  by  Hayya.  Observe 
that  these  values  range  from  .05  to  1.07,  with  a  median  value 
of  .36.  If  we  set  c  ■  .36,  then  (12)  yields  the  estimate  a  » 
7.7.  However,  evaluation  of  (8)  is  significantly  simplified 
if  a  is  integer,  for  then  the  Gamma  function  Y(a)  ■  (a-1)! 


'  1-10 

For  a«7,  c-1  /  /T  -.378,  while  a*8  gives  c«l  /  /F  •  .3523.  Hence, 
e  Gene  distribution  with  a-i  hes  e  coefficient  of  variation 
similar  to  the  median  c  value  observed  by  Hayya.  Now  suppose 
we  normalize  out  time  scale  so  that  E(t)  -  1.  Hence,  using  a 
■  8  in  (13)  yields  b  ■  1/8  -.IIS.  Finally,  substituting  these 
values  in  (8),  we  obtain 

(14)  g(t)  ■  _ 1 _  exp  (  -8t  ) 

V(8)  (1/8)® 

■  8  (8t)7  exp(  -8t) 

71 


which  yields 

(15)  g(t)  -  .0015873  (8t)7  exp(  -8t) 

Ne  have  used  this  equation  for  the  distribution  of  lead  times 
in  the  computer  code  presented  in  the  Appendix.  Finally,  suppose 
we  wish  to  estimate  the  unconditional  distribution  of  demand 
in  a  lead  time  using  the  P(  X  <  x  /  t)  distribution  defined  by 
(5),  (6)  and  the  lead  time  distribution  (15).  From  (7), 


(16)  P(  X  4  x  ) 


I 


x/R 


[o. 669  exp  (  .7979  z  )]g(t)  dt 


♦j  Cl.-  0.331  exp(  -  0.463  z  3  g(t)  dt 
'  x/R 


where  z  •  (x-Rt)/(QMAD  / 1  ). 


We  have  used  numeric  integration  to  evaluate  (16)  using  a 
step  size  of  dt  -  .1.  Our  computer  code  is  presented  in  the 
Appendix.  Subroutine  EXPLTD  (X,R,  QMAD,  QTRLT,  CUMPX)  computes 
the  cumulative  probability  CUMPX  *  P(  X  <  x)  for  given  values 
of  X,  R,  QMAD,  and  QTRLT,  where  QTRLT  denotes  the  expected  lead 
time  in  months,  and  the  other  terms  are  as  defined  above.  The 
MAIN  program  presented  in  the  Appendix  uses  subroutine  EXPLTD  to 
compute  the  unconditional  distribution  of  P(  X  <  x)  for  selected 
values  of  x  and  to  compute  and  print  associated  Laplace  distri¬ 
bution  values.  In  the  next  section,  we  present  the  results  of 
these  calculations. 


Section  II 


Sensitivity  Analysis 


To  compare  the  cumulative  distribution  functions  (CDF)  asso¬ 
ciated  with  the  Laplace  and  the  empirically-derived  model  of 
demand  in  a  leadtime,  we  developed  the  computer  code  shown  in 
Appendix  A*  We  then  used  this  code  to  evaluate  these  distribu¬ 
tions  for  a  number  of  parameter  sets.  Table  I I - 1  illustrates 
our  results  for  a  hypothical  item  with  a  demand  rate  R«300  units 
per  quarter,  a  demand  coefficient  of  variation  of  .2,  and  an 
average  leadtime  of  9  months.  In  the  table,  X  denotes  the  speci¬ 
fic  number  of  units  of  demand  in  a  given  lead  time,  while  the 
column  labeled  "EXPGAM"  shows  the  CDF  of  the  empirically  derived 
exponential -gamma  model.  That  is,  this  column  presents  the  cumu¬ 
lative  distribution  function  for  demand  in  a  lead  time  using 
the  exponential  model  for  forecast  errors  in  given  lead  time 
and  also  assuming  that  lead  times  are  independent  of  demand  and 
gamma  distributed  a  coefficient  of  variation  of  .353.  For  exam¬ 
ple,  comparing  these  two  columns,  we  see  that  there  is  an  801 
chance  that  demand  in  a  relenishment  leadtime  will  be  less  than 
or  equal  to  1330  units,  and  a  90.91  chance  that  demand  in  the 
leadtime  will  be  less  than  or  equal  to  1900  units. 

The  columns  labeled  "CONLT"  and  "LAPLACE"  represent  alternate 


cumulative  distribution  functions.  The  column  labeled  "CONLT" 


Where  X 

EXPGAM 

CONLT 

LAPLACE 

EDELT 


Table  II-l 

Cumulative  Probabilities 
for 

Three  Distributions  of  Lead  Time  Demand 


300.00 

CtW  * 

0.20 

LFAD  TIME  MONTHS 

*  9.00 

X 

EXPGAM 

CONLT 

LAPUCE 

EDELT 

EXF6-LAPL 

0. 

0.0004 

0.000! 

0.0000 

0. 

0.0003 

115.0 

0.0019 

0.0004 

0.0001 

0.1278 

0.0017 

230.0 

0.0093 

0.0011 

0.0004 

0.2556 

0.0099 

345.0 

0.0339 

0.0035 

0.0014 

0.3833 

0.0345 

460.0 

0.0972 

0.0099 

0.0046 

0.5111 

0.0926 

575.0 

0.1930 

0.0296 

0.0159 

0.6399 

0.1772 

690.0 

0.3195 

0.0692 

0.0536 

0.7867 

0.2650 

905.0 

0.4593 

0.2698 

0.1820 

0.8944 

0.2772 

920.0 

0.5965 

0.7039 

0.5958 

1.0222 

-0.0092 

1033.0 

0.8962 

0.0439 

0.9910 

1.1500 

-0.1849 

1150.0 

fi.TBSZ 

0.9177 

0.9650 

1.2779 

-0.178B 

1265.0 

0.954? 

0.9566 

0.9957 

1.4055 

-0.1349 

1390.0 

0.9020 

0.977i 

0.9970 

1.5333 

-0.0950 

1493.0 

0.9354 

C  9990 

0.9991 

1.651! 

-0.0637 

1610.0 

0.9394 

0.5937 

0  9997 

1.7699 

-0.0*14 

1725. ft 

0.9733 

0.9969 

! -9167 

-0.0295 

0.C3-? 

''.5993 

* 

7.  "M* 

-c.r.7! 

1935.0 

0  9931 

0.9991 

1.0000 

2.1722 

-0.0109 

2070.0 

0.9S35 

0.3533 

10000 

7.30C* 

-0.0069 

-  demand  in  the  leadtime  (units) 

-  Exponential 'Gamma  Model 

-  Exponential -Constant  Model 
•  Laplace  Model 

-  Standardized  lead  time  demand 


(Observed  lead  time  demand  X) 
(Expected  demand  rate)  (Expected  Leadtime) 


II 


is  the  emulative  distribution  function  for  deasnd  in  a  fixed 
lead  time  (i.e.,  there  is  no  variability  in  replenishment  lead- 
time),  using  the  exponential  approximation  to  the  distribution 
of  forecast  errors.  On  the  other  hand,  the  "Laplace"  column 
represents  the  corresponding  cumulative  probabilities  predicted 
by  the  LAPLACE  distribution  .  The  column  labeled  "EDELP"  is 
a  normalized  measure  for  demand  in  a  leadtime.  This  column  is 
obtained  by  dividing  the  number  of  units  X  demanded  in  a  leadtime 
by  the  expected  number  of  units  demanded  in  the  expected  lead- 
time.  In  this  case,  since  the  demand  rate  is  300  units  per  quar¬ 
ter  and  the  leadtime  is  9  months,  the  expected  demand  in  the 
expected  leadtime  is  900  units.  Consequently,  the  "EDELT"  column 
was  obtained  by  dividing  X  by  900.  Finally,  the  column  labeled 
"EXPG-LAPL"  presents  the  difference  between  the  culmulative  dis¬ 
tribution  functions  for  the  EXPGAM  and  LAPLACE  models. 

Figure  II-l  presents  a  plot  of  the  cumulative  distribution 
functions  in  Table  II-l.  In  the  figure,  the  solid  line  repre¬ 
sents  the  Laplace  cumulative  distribution  function,  while  the 
dashed  lines  present  the  constant  leadtime  and  exponential -gamma 
models,  respectively.  The  normalized  leadtime  demand  value  EDELT 
is  used  for  the  X-axis  in  this  plot.  As  shown  in  the  figure, 
there  are  significant  differences  between  the  Laplace  and  Expo¬ 
nential-Gamma  models.  For  example,  if  one  wishes  to  achieve 
a  90  percent  fill  rate,  the  EXPGAM  model  indicates  that  the  safety 
stock  should  be  set  to  1.53  times  EDELT,  the  expected  demand 
in  the  expected  leadtime,  or  1380  units.  On  the  other  hand. 


the  Laplace  model  indicates  that  the  901  fill  rate  may  be  achieved 
with  a  safety  stock  of  about  1.2  tiaes  BDELT,  or  1080  units. 

On  the  other  hand,  the  two  curves  cross  at  approxiaately  the 
EDBLT*1  value,  and  the  cuaulative  distribution  functions  of  all 
three  functions  are  very  siailar  in  this  region. 

Figure  II-2  presents  a  plot  of  the  three  distributions  for 
a  case  in  which  deaand  per  quarter  is  300  units  and  leadtiae 
is  9  aonths,  but  the  coefficient  of  variation  of  deaand  per  quar¬ 
ter  has  been  increased  to  .5.  Notice  that  the  three  curves  are 
closer  in  this  case,  but  that  there  are  still  substantial  differ¬ 
ences  aaong  the  curves.  Figure  I I - 3  presents  a  siailar  plot 
when  the  coefficient  of  variation  of  demand  per  quarter  is  .8. 

The  curves  are  now  even  more  similar  than  for  the  .5  case,  but 
significant  differences  aaong  the  curves  still  exist,  particu¬ 
larly  in  the  801  and  above  fill  rate  region. 

To  obtain  further  insights  into  the  relative  behavior  of 
these  three  curves,  we  ploted  a  number  of  other  combinations 
of  parameters.  In  our  first  sensitivity  study,  we  were  inter¬ 
ested  in  the  effects  of  changes  in  item  demand  rate  and  demand 
variability  upon  the  overall  shapes  of  the  curves.  Our  results 
are  presented  in  Figures  I I - 1  through  11-12,  while  the  specific 
parameter  sets  investigated  are  shown  in  Table  1 1 - 2 .  As  shown 
in  the  Table,  we  developed  curves  for  demand  rates  of  300,  30, 

3,  and  .3  units  per  quarter,  respectivily,  and  for  coefficients 
of  variation  of  demand  per  quarter  of  .2,  .5,  and  .8.  In  all 


SISTRtftUTION  FUNCTION* 


Table  II -2 

Parameter  Sets  for  Deaand  Rate  and  Deaand  Variability  Sensitivity  Tests 

Coefficient  of  Average 


Figure 

Deaand 

per  Qtr 

Variation  of 
Deaand  per  Qtr 

Replenishment 

Leadtiae  (Months) 

II-l 

300 

.2 

9 

II-2 

IV 

.5 

9 

II-3 

VI 

.8 

9 

II-4 

30 

.2 

9 

II-5 

ft 

.5 

9 

II-6 

If 

.8 

9 

II-7 

3 

.2 

9 

II-8 

VV 

.5 

9 

II-9 

.8 

9 

11-10 

.3 

.2 

9 

11-11 

•1 

.5 

9 

11-12 

ft 

.8 

9 

of  these  cases,  we  assumed  that  replenishment  leadtime  was  gamma 
distributed  with  a  coefficient  of  variation  of  .353  and  an  aver¬ 
age  leadtime  value  of  9  months. 


As  may  be  seen  in  Figures  I I- 4  through  11-12,  as  the  coef¬ 
ficient  of  variation  of  demand  per  quarter  increases,  differences 
among  the  three  leadtime  demand  curves  diminish.  The  greatest 
differences  among  the  curves  occur  at  low  values  of  the  coeffi¬ 
cient  variation,  and  the  differences  decrease  as  the  coefficient 
of  variation  increases.  However,  even  for  coefficient  of  varia¬ 
tion  values  of  .8,  significant  differences  among  the  curves  exist 
for  fill  rates  in  the  801  or  higher  range. 


Sensitivity  to  Average  Replenishment  Leadtime 

We  also  developed  a  number  of  plots  to  investigate  the  sensi¬ 
tivity  of  the  three  leadtime  demand  curves  to  changes  in  the 
average  replenishment  leadtime.  Table  II-3  summarizes  the  para¬ 
meters  sets  used  while  Figures  11-13  through  11-24  present  our 
results.  In  this  case,  we  investigated  average  replenishment 
leadtimes  of  4,  9,  and  12  months,  while  demand  per  quarter  was 
set  to  300,  30,  3  and  .3  units  per  quarter,  repectively.  In 
all  of  these  calculations,  the  coefficient  of  variation  of  demand 
per  quarter  was  sent  to  .5  while  the  coefficent  of  variation 
of  replenishment  leadtime  was  set  to  .353. 


Table  II -3 


Parameter  Sets  for  Replenishaent  Leadtiae  Sensitivity  Tests 


Deaand 
Figure  per  Qtr 


IT- IS 

— 5TO 

11-14 

ft 

11-15 

«• 

11-16 

30 

11-17 

ft 

11-18 

ft 

11-19 

3 

11-20 

tt 

11-21 

ft 

11-22 

.3 

11-23 

tt 

11-24 

tt 

Coefficient  of 
Variation  of 
Deaand  per  Qtr 

75 - 


.5 


.5 


.5 


Average 
Replenishaent 
Leadtiae  (Months) 

75 - 

9 

12 

6 

9 

12 

6 

9 

12 

6 

9 

12 


Figures  11-13  through  11-24  present  relationships  aaong  the 
curves  which  are  very  siailar  to  those  observed  in  Figures  II- 
1  thru  11-12.  As  the  leadtiae  increases,  slight  changes  in  the 
curves  take  place,  but  these  are  hard  to  observe  in  the  graphs. 
In  all  cases,  significant  differences  exist  aaong  the  curves 
in  the  80+1  fill  rate  ranges. 
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Appendix  a 

Progress  for  Laplace 
and 


Bapirical  Approximation  Calculations 


0 


1 


10*EXPLTD.S— COMPUTE  P(X«<  x>  FOR  EXP.  APPROX.  ANI  GANNA  LEADTINE 
2#C 

30C  ASSURE  1  PERIOD  *  E1PECTEI  LEADTIME 
ARC  LEADTINE  HAS  MEAN  ■  1. 

SOC  COEFFICIEIT  OF  VARIATION  -  .153 

AOC 

7IC  »  »  EXPECTED  IE I AND  IN  LEADTINE 

SIC  816  >  STAIDARD  DEVIATION  OF  DEMAND  IN  EXPECTED  LEADTINE 

PIC  RNAD  >  NAD  OF  DEMAND  IN  EXPECTED  LEADTINE 

1IOC 

110C  CUNPHCUN.  PROB  OF  LEADTINE  DISTRIBUTION 

120C  CUNPX*  CUN.  PROD  OF  DEIAND  X  IN  LEADTINE 

1  SOC 

140C  DT  >  TINE  DELTA 

190C  IT  >  P  (LEADTINE  *  T) 

1AOC  ZT  *  STANDARDIZED  ERROR 

HOC 

100  CONNOR/ IUT / I NT ( 21 ) 

HOC 

200C  SET  PRINT  FLASS 
210C 

221  I NT ( 1 0  > ■ I BETL 

230  IITM1HIPNTSZ 

240  CALL  FPARAN< 1,132) 

2S0C 

260  PRINT, "OUTPUT  TO  FILE  08*  < 1*101" 

270  READ.IOUT 

20OC 

2V0C 

30IC 

310  PRINT, "PRINT  DETAILS  T  DETAIL  STEP  SIZE*" 

320  READ  ,  1DETL,1PNTSZ 

330  IFIIPNTSZ.LT. 1 )  IPNTSZ*1000 

340C 

342  S  CONTINOE 

3S0C  SET  PRINT  FLAGS 

3A0C 

371  INT(1l)»IDETl 

300  IIT(1 1 )  * I PNTSZ 

390  PRINT, •REA#  AND  COFV  IF  DEMAND  PER  OTR,  AND  IT  IN  MONTHS" 

410  READ, R, COFV,  RLT 

405  RNAD  *  0.8*C0FV*R 

410  PRINT  23, R, COFV,  RLT 

420  23  FORMAT <//*R  ■■,F8.2,a  COFV  «",F|.2, 

430D  "  LEAD  TINE  NONTNS  **,F8.2//,'i 

440C 

450  IF( IDETL.6T.0) 

4001  PRINT,"  X  I  ST  PXf  CVMPT  LVNPX’. 

4701  "  ZT  PZ" 

4  IOC 


EXPLTD 


CONFUTE  LEAST 1NE  II  IUARTEIS 
ESTIMATE  NAD 


4  TOC 
SOOC 
510C 
S20C 
S30C 

541  QTRLT*RLT/3. 

5 SO  TQTRaQTILT 

SSS  EILT  »TITI*R 

5  SAC 

540C  EVALUATE  P<  X  <*  X  )  FOI  1  *0  TO  MEAN  *  1*KMND 

570C 

510  XNAX>R«0TRLT+1. •RIAD*SQRT< STILT ) 

5V0C 

400  OX  «  INAX/1 I. 

410  IFIIX.LT.  1.)  IX  >1. 

420  1F10X.ST.  1.)  OX  *IFIX<IX  ♦  0.5) 

430C 

440C 

450C  INITIALIZE  VARIABLES 

440C 

670C 

480C  CONFUTE  CUHPX«P( X<*X> 

6  TOC 

700C - - - - - 

7 tOC  SCSI I  "X"  EVALUATION  LOOP 

720C - - 

730  PRINT,"  X  EXF6AN  CONLT  LAFLACE  EIELT* 

7351  EXP6-LAPL" 

740  PRINT,"  " 

750C 

740  X*0. 

770  00  200  IX* 1,1 00 

7I0C 

790  CALL  EXPLTD(X(R,RIIAI,aTlLT,CUnPX> 

800C 

810C  COMPOTE  FIXES  LEADTIME  PROBABILITY 

820C 

830  Z><X-l*TQTR)/<RMAD*SORT<IQTR)) 

840  IF<Z.LE.O)CPFLT*.4A9*EXP<0.7?7?«Z1 

830  IF<2.6T.0.)CPFLT*1.-.331*EXP<-.4A3»Z> 

8A0C 

870C  COIPOTE  LAPLACE  PNQBABTLITT 

8S0C 

80S  SI 8» .5V45*RNAI* ( . 1235*0. 42*25*AL T ) 

8 VO  RK-(X-R*TQTR>/<  SIGI 

900  IF(»K.LE.O.>CPLPC*.5*EXP<1.«142>*RK) 

910  IF(RK.6T.0.)CPLPC*1.-.5=*EXP(  -1.4147*RX> 

920C 

930C  PRINT  RESULTS 

940C 

942  XN  •  X/EDLI 

945  DIFF»CUNPI-CPtPC 

950  PAINT  123,X,CU«FX,CPFLT,CPLPC,XN,I1FF 

9A0  123  FORMAT  <F 12.1, 5F 10.4) 


EXPLTD 


970C 

980  LINE-LIME  ♦  1 

990  IF(I0IT.GT.0>URITE<8,133>LINEfX,CUMPXrCPFlT,CPLPC,XN 

1900  133  FORMAT < 15, F8. 1 ,4F1 0.4) 

1010C 

1030C 

1030C  INCREMENT  X  AMD  CHECK  IF  DONE 

1040C 

1050  X*X+DX 

1060  IFICNMPX.MT .  0.99)  SO  TO  300 

1070C 

1000C 

1090C - ENI  OF  "X*  LOOP 

1100  200  CONTINME 
1  HOC 

1120  300  CONTINUE 

1121  LINE  «  L1NEM 

1122  IF(IOUT.ST.O)  UNITE (8,1 J3)IIME,0.,0.,0.,0. ,0. 

1123C 

1130C 

1140  MRITE( 4,423) 

1145  423  FORMAT! /////“CONTINUE  »  <1*TES)"/> 

1150  REAM, ICON! 

1140  I F ( ICOMT.NE.  1)  STOP 

1170  60  TO  5 

1 180  END 

1190  SU0R0UT1NE  EXPLTDU  ,RfartADfQTRlT  .CUMPK  > 

1200C 

1210  CONH0N/IMT/IMT(2O) 

1220C 

1230C  IET  WRITE  FLAMS 

1240C 

1250  IDETL  *  IMT ( 1 0 ) 

1240  IPNTS2  *  IUT<in 

1270  CUMPT*0. 

1280  CUHPX=0. 

1290C 

1300C  INITIALIZE  PDF  PARAMETERS 

1 31 OC 

1320  A1 *  0.331 

1330  II  *  -0.443 

1340C 

1350  A2*  0.649 

1340  12 «  0.7979 

1370C 

1 JOOC  SANAA  CONS  I  AN  r  FOR  AEAN  -1  ANO  CV=.I5< 

1390  C1«  0.0015873 

1400C 
1 41  OC 


HXPLTD 


142IC - - - - 

1430C 

1440C 

145IC  INITIALIZE  FOR  T  INTEGRATION 

1 440C 

1470  DT  -  .1 

1410  I  >  IT 

1490  CUNPT*!. 

ISO*  CUNPX-4. 

15I0C 

1520C  IESIN  "T"  INTEGRATION  10NP 

1S30C 

1 S40  10  101  I«1 r  tOO 

1550C 

1540C  CONFUTE  STANDARI1ZED  ERRNR  ZT 

1370C 

1580  TOTR*  TtQTRLT 

1590  IT  ■  («  -  R*TQTR)/  (QNAMIQRT ( TQTR)  ) 

1400C 

141 OC  CONFUTE  p(  Tl 

1420C 

1434  ST  «  C1*(I.*T >*•?•  EXP(-I.*T)  •  IT 

1440C 

1450C  CINPNTE  9(1  <•  ZT  J  T» 

1440C 

147#  IF*  ZT.LE.O.)  PZ  «A2*EXP(B2*ZT> 

1480  IF (  ZT.GT.4.)  PZ*  1.  -  4HEIPI  B1*ZT) 

1490C 

1700C  CONFUTE  P<  X  <«  i  !T>  P<  I) 

1 71  OC 

1720  PIT  -  PZ«GT 

1730C 

1740C  UPDATE  CUMJIATIVE  PROBABILITIES 

17S0C 

1740  CUNPT  *  CUNPT  ♦  GT 

1770  CUNFX  =  CUNPX  ♦  FXT 

1 780C 

1790  IPRNTM 

1800  IF<NOD(I,IPNTSZ).EQ.O>  IPRNT*1 

1110  IF< IBETL.LE.O)  IFRNT>0 

1120  IF(IPRNT.ST.O)  UR1TEU,4J>  X,T ,GT , FXT, CUNFT, CUNFX fZT,9Z 

1130  43  FORMAT ( 2F8 . 2, 4F 1 0 . 4 ) 

1I40C 

1850C 

1840C  INCREMENT  T 

1I70C 

1810  T  »  T+IT 

1890C 

19O0C  STOP  IF  CUNPT  >  .991 

19I0C 

1920  IF (CUNPT .ST .  0.999)  GO  TO  120 

1930C 

1940C - ENI  OF  "T"  I  OOP - 

1950  100  CONTINUE 

1940C 

1970  120  CONTINUE 
1980  RETURN 


1990  .  ENI 


EXPLTD 


